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The well-known results concerning a dilute Bose gas with the short-range repulsive interaction 
should be reconsidered due to a thermodynamic inconsistency of the method being basic to much 
of the present understanding of this subject. The aim of our paper is to propose a new way of 
treating the dilute Bose gas with an arbitrary strong interaction. Using the reduced density matrix 
of the second order and a variational procedure, this way allows us to escape the inconsistency 
mentioned and operate with singular potentials of the Lennard-Jones type. The derived expansion 
of the condensate depletion in powers of the boson density n = N/V reproduces the familiar result, 
while the expansion for the mean energy per particle is of the new form: e = 2irh 2 an/m{l + 
128/ (15-s/tt)"^ na 3 (1 — 56/8a) + •■•}, where a is the scattering length and b > stands for one more 
characteristic length depending on the shape of the interaction potential (in particular, for the hard 
spheres a = b). All the consideration concerns the zero temperature. 
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It is well-known that to investigate a dilute Bose gas of 
particles with an arbitrary strong repulsion (the strong- 
coupling regime), one should go beyond the Bogoliubov 
approach (weak-coupling case) and treat the short- 
range boson correlations in a more accurate way. An 
ordinary manner of doing so is the use of the Bogoliubov 
model with the "dressed" , or effective, interaction poten- 
tial containing "information" on the short-range boson 
correlations (see Ref. ||). Below it is demonstrated that 
this manner leads to a loss of the thermodynamic consis- 
tency. To overcome this trouble, we propose a new way of 
investigating the strong-coupling regime which concerns 
the reduced density matrix of the second order (the 2- 
matrix) and is based on the variational method. 

The 2-matrix for the many-body system of spinless 
bosons can be represented as j|: p2(r'i, r 2 ; r 2 ) = 
-FMri, r 2 ; r[, r' 2 )/{N(N — 1)}, where the pair correlation 
function is given by 

F 2 (r 1 ,r 2 ;r' i y 2 ) = (^ (v^ ^{v'^v',)) . (1) 

Here 4>( r ) and V^( r ) denote the boson field operators. 
Recently it has been found j|,[3| that for the uniform sys- 
tem with a small depletion of the zero-momentum state 
the correlation function (|l|) can be written in the ther- 
modynamic limit as follows @||: 

^ 2 (ri,r 2 ;ri,r 2 ) = n 2 tp*(r)tp(r') 

+ 2n ° j (Ss Hq ^q/2( r Vq/2(0 exp{*q(R' - R)}, (2) 

where r = ri — r 2 , R = (ri + r 2 )/2 and similar rela- 
tions take place for r' and R', respectively. In Eq. (^) 
n = Nq/V is the density of the particles in the zero- 
momentum state, n q = (a q a q ) stands for the distribu- 
tion of the noncondensed bosons over momenta. Besides, 
ip(r) is the wave function of a pair of particles being both 
condensed. In turn, ipq/2 (r) denotes the wave function 



of the relative motion in a pair of bosons with the total 
momentum ?iq, this pair including one condensed and 
one noncondensed particles. So, Eq. (||) takes into ac- 
count the condensate-condensate and supracondensate- 
condensate pair states and is related to the situation of a 
small depletion of the zero-momentum one-boson state. 
For the wave functions <p{r) and f P (j) we have 

<p(r) = 1 + if>(r), ¥>p(r) = y/2 cos(pr) + ^(r) [p + 0) (3) 

with the boundary conditions ip(r) — > and ipp( r ) ~ * 
for r — * oo. The functions ip(r) and ^p( r ) can explicitly 
be expressed in terms of the Bose operators ap and a p 
[0]. In particular, 

i>{k) = J tp(r) exp(-ikr) d 3 r = (aw a_k)/no • (4) 

Having in our disposal the distribution function rik and 
the set of the pair wave functions <p(r) and c/? p (r), we are 
able to calculate the main thermodynamic quantities of 
the system of interest. In particular, the mean energy 
per particle is expressed in terms of nj~ and g(r) via the 
well-known formula 

where Tk — h 2 k 2 /2m is the one-particle kinetic energy, 
n = N/V stands for the boson density and the relation 

g(r) = F 2 (ri,r 2 ;ri,r 2 )/n 2 . (6) 

is valid for the pair distribution function g(r). 

The starting point of our investigation is the weak- 
coupling regime which implies weak spatial correlations 
of particles and, thus, is characterized by the set of the 
inequalities 



1 



|V(r)|«l, |V> P (r)|«l 



(7) 



Specifically, the Bogoliubov model corresponds to the 
choice HI 



|V(r)|«l, ^ P (r)=0 



(8) 



Besides, owing to a small depletion of the Bose conden- 
sate (n — no)/n we have for the one-particle density ma- 
trixFi(r) = (V> f (ri)V(r2)): 



Fi(r) 



d 3 k 



n k 



(2tt) 3 n 



exp(ikr) 



n — nn 
< < 1. 



So, investigating the Bose gas within the Bogoliubov 
scheme, we have two small quantities: ip(r) and F\{r)/n. 
This enables us to write Eq. (||) with the help of (||) as 
follows: 



g(r) = 1 + 2tp(r) + 



d 3 k 
n J (2tt) 3 



n k exp(ikr) 



(9) 



where we restricted ourselves to the terms linear in tp(r) 
and Fi(r)/n and put ip*{r) — 4 , ( r ) because the pair wave 
functions can be chosen as real quantities. Equations for 
tp(k) and n k can be found varying the mean energy (||) 
with (^) taken into account. However, previously one 
should realize an important point, namely: n k and ip(k) 
can not be independent variables. Indeed, when there 
is no interaction between particles, there are no spatial 
particle correlations either. So, ip(k) — and, since the 
zero-temperature case is considered, all the bosons are 
condensed, n k = 0. While "switching on" the interaction 
results in appearing the spatial correlations and conden- 
sate depletion: ^(fc) ^ together with n k ^ 0. In the 
framework of the Bogoliubov scheme ip(k) is related to 
n k by the expression 



n k (n k + 1) = nlip 2 (k) 



(10) 



Indeed, the canonical Bogoliubov transformation JlJ im- 
plies that 

flk = u k a k + UfcaLk, al = u k al + v k a-k, (11) 
where 



2 2i 

ut—vt = 1. 



(12) 



At zero temperature (a^ak) = and, using Eqs. (Q) and 
(O) we arrive at 



(13) 



n k = vl, tp(k) = u k v k /n a . 



With Eqs. ([12]) and (]13j) one can readily obtain ([UJ). 

Now, let us show that all the results on the thermody- 
namics of a weak-coupling Bose gas can be derived for 
the Bogoliubov scheme with variation of the mean energy 
(||) under the conditions (||) and (|l0|). Inserting @ into 



and, then, varying the obtained expression, we arrive 



at 
Se 



d 3 k 



; T k + n$(fc) ) — + n®(k)5$(k) 



(T k + n$(fc)) 



(14) 



Relation ( |l0| ) connecting 4>(k) with n k results in 
(2n k + l)6n k $(k) f d 3 q 



5ip(k) 

2nl^{k) 

where the equality 

n = n + 



"o 



(2ny 



d 3 k 

(>j3 



n k 



Sn q , (15) 



(16) 



is taken into consideration. Setting Se = and using 
Eqs. (|TJ) and @, we derive the following expression: 



2T fc V(fc) = — *(fc)(l + 2n fc ) 



I — 77 

+2nM) $(fc) + — 
\ n 



d 3 q 
(27T) 



Here one should realize that Eq. (17) is able to yield re- 



sults being accurate only to the leading order in (n — 
no)/n because the used expression for ofr) given by (|9 ) 
is valid to the next-to-leading order f|. So, Eq. (J17 ) 
should be rewritten as 



-2T k ip(k) = $(fc)(l + 2n k ) + 2ni>(k)<S>(k). (18) 

Eq. ( 1^) is an equation of the Bethe-Goldstone type or, 
in other words, the in-medium Schrodinger equation for 
the pair wave function. As 2^(k)(n k +mp(k)) is the prod- 
uct of the Fourier transforms of $(r) and n(g(r) — 1), we 
can rewrite Eq. (jig) in the more customary form 



m 



VVW = «(r) + n J $(|r - y|) (g{y) - l) d 3 y. (19) 

The structure of Eq. (|I7J) is discussed in the papers HQ . 
Here we only remark that the right-hand side (r.h.s.) of 
( |l9| ) is the in-medium potential of the boson-boson inter- 
action in the weak-coupling approximation. The system 
of equations (JTJJ) and ( |l8|) can easily be solved, which 
leads to the familiar results [Q: 



(20) 




2\/T? + 2nT k $(k) 



Now we are able to demonstrate that the investiga- 
tion of the strong-coupling case based on the Bogoliubov 
model with the effective boson-boson interaction, results 
in a loss of the thermodynamic consistency. Indeed, as 
it was shown in the previous paragraph, any calculating 



2 



scheme using the basic relations of the Bogoliubov model 
@j (0) conclusively leads to Eqs. (|i^)-(|20|) provided this 
scheme does yield the minimum of the mean energy. In 
this case Eqs. (18)-(P0[) certainly includes the quantity 
$(r) which is the "bare" interaction potential appear- 
ing in (g). The use of the Bogoliubov model with the 
effective interaction potential substituted for $(r) can 
in no way disturb the relations given by (g) and (10). 
And Eq. (g) is the same in both the weak- and strong- 
coupling regimes. Thus, any attempts of replacing <I>(r) 
by the effective "dressed" potential without modifications 
of (g) and © results in a calculating procedure which 
does not really provide the minimum of the mean energy. 
It is nothing else but a loss of the thermodynamic con- 
sistency. Remark that we do not mean, of course, that 
the t-matrix approach or the pseudopotential method 
can not be applied in the quantum scattering problem. 
It is only stated that the usual way of combining the 
ladder diagrams with the random phase approximation 
faces the trouble mentioned above. Though our present 
investigation is limited by the consideration of the Bose 
systems, the derived result gives a hint that the similar 
situation is likely to take place in the Fermi case, too. In 
this connection it is worth noting the problem associated 
with the lack of self-consistency of the standard method 
of treating the dilute Fermi gas M . 

The strong-coupling regime is characterized by sig- 
nificant spatial correlations. So, Eq. (g) resulting in 
(g) is not relevant for an arbitrary strong repulsion be- 
tween bosons at small separations when we have ?/>(0) = 
-1, V P (0) = -V2 (see Refs. |,|). Therefore, to in- 
vestigate the strong-coupling regime, Eq. (g) should be 
abandoned in favor of (g). Expression (g) is accurate to 
the next-to-leading order in (n — no)/n. So, using (g) 
and (o), we can write 



g(r) =<p (r) + - 
n 



(27T) 
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(<4 /2 (r)-^ 2 (r)). ( 



21) 



Let us now perturb t[>(k) and n(k). Working to the first 
order in the perturbation and keeping in mind conditions 
© and d|), from (g) we derive: 



2T k ^{k) = U(k)(l + 2n k ) + 2n^(k)U'(k) (22) 



with 



and 



U(k)= tp(r)®(r)exp(-ik.r)d 3 r 



U'{k)= / Ul /2 (r)~v 2 (r) $(r)d 3 r 



(23) 



(24) 



Using Eqs. ([2 3D, (|24j) as well as the relation V'k(r) — > 
y/2ip(r) (k —> 0) (see the boundary conditions (||)), 
we obtain U(0) ^ U'(Q). This implies that the system 
of Eqs. (|io|) and (||5) is not able to yield the relation 



n k oc 1/fc (fc — > 0) following from the "1/fc 2 " theorem 
of Bogoliubov for the zero temperature || . Indeed, let 
us assume n k — > oo for fc — > 0. Then, from Eq. (10) at 
n = rig we find n\ip(k)\/n k — > 1 when k — > 0. On the 
contrary, Eq. @ gives n${k)\/n k -> U(0)/U'(Q) ^ 1 
for k — > 0. So, consideration of the Bose gas based on 
Eqs. (g) and ( |Io| ) does not produce satisfactory results. 
Nevertheless, it is worth noting that Eq. (|2^) has an im- 
portant peculiarity which differentiate it from Eq. jl^ ) 
in an advantageous way. The point is that in both the 
limits n — > and k — * oo Eq. (|22| ) is reduced to 



— V 2 (p(r) + $(r)ip(r) = 0. 



(25) 



As it is seen, this is the exact "bare" (not in-medium) 
Schrodinger equation, other than its Born approximation 
following from ([l9]). Thus, we can expect the line of our 
investigation to be right. 

As it was shown in the previous paragraph, an ap- 
proach adequate for a dilute Bose gas with an arbitrary 
strong interaction can not be constructed without mod- 
ifications of Eq. ([l0|). This is also in agreement with a 
consequence of the relation 



|(a k a_ k )| 2 < (a k aj c )(aL k a_ 



(26) 



resulting from the inequality of Cauchy-Schwarz- 
Bogoliubov § |(AB)| 2 < {AA^){B^B). With (g) and 
( p6| ) one can easily derive ng^ 2 (fc) < n k {nk + 1). Thus, 
it is reasonable to assume that Eq. ([l(]) takes into ac- 
count only the condensate-condensate channel and ig- 
nores the supracondensate-condensate ones. Now the 
question arises how to find corrections to the r.h.s. of 
Eq. ([l0]) . At present we have no regular procedure allow- 
ing us to do this in any order of (n — no)/n. However, 
there exists an argument which makes it possible to re- 
alize the first step in this direction. The matter is that 
the alterations needed have to produce the equation for 
?/>p(k) which is reduced to the equation for ip(k) in the 
limit p — > 0. Though this requirement does not uniquely 
determine the corrections to Eq. (|Io|), it turns out to 
be significantly restrictive. In par ticular, even the sim- 
plest variant of correcting Eq. (|l(j| ) in this way, leads to 
promising results. Indeed, this variant is specified by the 
expression 



n k (n k + 1) = nl V )2 (fc) + 2n 



d 3 q 
(2tt) 3 



n,^ 2 (k). (27) 



Eq. (|27]) is valid to the next-to-leading order in (n — 
no)/n. So, we may rewrite it as 

n k (n k + 1) = n 2 v7 2 (fc) + 2n [ ^n 9 (^ 2 /2 (k) - ^(k)) . 



(28) 
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Perturbin g il> (k) and n k and bearing in mind conditions 
( pl| ) and (||§), (^) gives Eq. (|2^) again. However, now 
C/'(fc) obeys the new relation 



U'(k) = 



r (yl /2 {v)-V 2 {r))Hr)d?r 

I- d 3 q U(q)^l /2 (cQ-^(q)) 
J ( 27r ) 3 i>(q) 



(29) 



which significantly differs from (24j). Indeed, the choice 
of the pair wave functions as real quantities implies that 
operating with integrands in (^3) and ( |29]) , one can ex- 
ploit - 0p( r ) — V2ip(r) oc p 2 at small p [^0|. For k — > 
this provides £/'(fc) — J7(fc) = £fc = c/c 4 + • • • . Similar 
to Eq. (|l8|), Eq. p2] ) can yields results correct only to 
the leading order in (n — no)/n. So, it has to be solved 
together with ([l0|) where n§ should be replaced by n 2 , 
rather than with J2£|). This leads to the following rela- 
tion: 



1 



Tfc + nU(k) 
If* + 2nf k U(k) 
U(k) 



1 



(30) 
(31) 



2a/T 2 + 2nT k U{k) 



where T k = T k + nt k ■ In the limit k — > Eq. ([H]) gives 

~ (■\JnmU(0)/hk — l)/2, which is fully consistent 
with the "1/fc 2 " theorem of Bogoliubov for the zero tem- 
perature ||. Eqs. ( p3| ) and ( pi] ) should be solved in a 
self-consistent manner. So, for n — > one can derive 



U(k) = U {Q) (k)(l + 8VnaVV5r). 



(32) 



Here U^°\k) — J v?^(r)$(r) exp(— ik.r)d 3 r, where 
^j(°)(r) obeys Eq. (|25|). Further, substituting fc = \fny 
in the integral for the condensate depletion (n — na)/n = 
1/(2tt) 3 / c , dfc4^fc 2 n fc /n, we obtain the familiar result 



(n — no)/n = 8Vna 3 / + 



(33) 



a being the scattering length. Inserting (pl|), ( |30| ) and 
( |3l| ) into Eq. (|^) and using (^), in a similar manner we 
derive 



L50F 



(34) 



where b > is one more characteristic length defined as 



h / (V ^ 1 " ( '- i)J ' / ' : '-- 



(35) 



As it is seen, the well-known result of papers j| can 
be derived from ([34]) with the choice 6 = 0. How- 
ever, this approximation is rather crude because the case 
of the hard-sphere interaction ($(r) = (r > a) and 
<I>(r) — > oo (r < a)) is specified by 6 = a: 




(36) 



In the general case and, in particular, for the singular 
potentials of the Lennard-Jones type we have a ^ b. Re- 
mark that the last term in the r.h.s. of ( |27| ) does not 
make any contribution into the results given by (|33| ) and 
(|34|). However, the next orders in the expansions of the 
energy and depletion depend on its contribution essen- 
tially. 

Concluding let us take notice of the important points 
of this Letter once more. It was demonstrated that ther- 
modynamically consistent calculations based on (j^ and 
(|Io| ) conclusively result in Eqs. (|l8|)-(po|). Therefore, us- 
ing the Bogoliubov model with the "dressed" interaction 
does not provide the satisfactory solution of the prob- 
lem of the strong-coupling Bose gas. As it was shown, 
when investigating this subject, one should go beyond 
the Bogoliubov scheme. To do this, we developed the 
approach reduced to the system of Eqs. (p3|), (£9|), ( |30| ) 
and (pl|) . This equations reproduce the well-known re- 
sult ([33|) for the condensate depletion and yields the new 
expansion (34|) in powers of n for the energy, ( |36] ) being 
the particular case of the hard spheres. One can expect 
alterations for the excitation spectrum, too. 
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